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SUMMARY 


A numerical calculation has been carried out to evaluate the 3x3 
cross section matrix involved in the electron impact excitation of the 
ground state of H atom to the Zs and the 2p levels. The method of solu- 
tion is that of atomic eigenstates expansion. Similar calculation has been 
carried out by Burke, Schey and Smith. In this paper, however, the 
definite integral terms in the coupled radial differential equations are 
eliminated through some linear transformation of the radial functions , 
thus avoiding iteration of these equations. As a result, the equation of 
reciprocity and the equation of continuity of currents are numerically 
satisfied with an error to value ratio less than 1 per 1000 on the aver- 
age, and the maximum of this ratio, except for a few cases, has been 
kept below 1%. A simple perturbation theory has been developed to 
evaluate the effect of the long range and the centrifugal potentials. The 
five cross sections, ls-Zs, ls-Zp, ls-ls, Zs-Zs and Zp-Zp,are tabulated. 

The agreement of the ls-Zs cross sections with the experimental 
results are satisfactory while the calculated ls-Zp cross sections are 
higher than the corresponding experimental values. The Zs-Zs and the 
2p-Zp cross sections have large values, specially at the thresholds. A 
Ramsauer effect in the partial cross sections of the Zs-Zs elastic 
scattering is observed. By comparison with the eigenstate expansion 
calculation it isfoundthat the Born approximation, despite its simplic- 
ity, gives meaningful results for low and close to threshold energies of 
the bombarding electrons. In this paper the effect of the electron spin 
and the exchange potentials are also investigated. 
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I. INTRODUCTION 


The recent experimental results of Lichten and Schultz 1 on one hand, 
and Stebbing et al. 2 on the other, of the Is — ► 2s transition cross section 
in the hydrogen atom by electron impact and the apparent discrepancy in 
these measurements, and the measurement of Is — *2s transition by Fite 
et al. 3 also in hydrogen, necessitate accurate calculation for these cross 
sections. 

Calculation of the excitation cross sections in atomic hydrogen by 
electron impact corresponds to the solution of the problem of three in- 
teracting bodies, one proton and two electrons. By taking the position 
of the proton as the center of mass, the problem will reduce to the task 
of finding the nonseparable wave function of the system of the two elec- 
trons with an attractive center of force. Such solution has not been 
found yet. However, if this wave function is expanded in terms of the 
eigenstates of the hydrogen atom, it is shown in this paper that through 
numerical integration the coefficients of the expansion, which are func- 
tions of the position vector of the free electron, can be found exactly. 
When an infinite number of terms are included in the expansion the solu- 
tion to the problem is exact. Furthermore, the expansion has the ad- 
vantage that the asymptotic form of its coefficients are automatically the 
asymptotic form of the free electron wave function scattered from dif- 
ferent atomic states, which are simply related to the excitation cross 
sections. 
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In this paper atomic states Is , 2s, 2p are included in the expansion, 
and by antisymmetrizing the two electron wave functions according to the 
exclusion principle, some contribution from the continuum in the expan- 
sion is also taken into account. The first calculation of this type was 
performed by Marriot 4 whose expansion consisted of the Is and the 2s 
states in order to calculate the Is — *2s transition cross section. This 
calculation was extended by Smith 5 to higher total orbital angular mo- 
menta of the system. Percival and Seaton 6 have formulated the eigen 
state expansion technique in general, and have tabulated the coefficients 
of the integro-diff erential equations for s, p and d atomic electrons. 
vVhile this paper was in preparation, Burke, Smith and She y, 7 * 32 using the 
tables of reference (6) for three states Is, 2s, 2p, have integrated the 
resulting integro-diff erential equations. Their technique of solution 
is substantially different from that of this paper in more than one re- 
spect, and it is believed that the results presented here are more 
accurate. 

The degeneracy of the 2s and the 2p levels of the hydrogen atom 
makes it necessary, as pointed out by Massey, 8 that any calculation 

concerning Is — >2s or Is >2p transitions contain these two states 

simultaneously. In this respect the present calculation and that of refer- 
ence (7) is superior to that of reference (4), where Is and 2s states are 

included to calculate Is >2s transition, or that of Khashaba and Massey, 9 

where Is and 2p states are included to calculate Is »2p transition. 
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We treat the present problem in the following manner. In Sec. IIA 
we derive the integro-differential equations of the problem independent 
of reference (6). Sec. IIB is the derivation of the transmission matrix 
and the cross sections. In Sec. IIC we derive a useful relationship for 
numerical integration based on the symmetry of the interacting potentials. 
Sec. IID is the derivation of the partial wave Born approximation trans- 
mission matrix. Different parts in Sec. Ill deal with techniques of 
numerical integration and the effects of the long range potentials. The 
transformation of the radial wave functions presented here avoid the 
need for iteration of the coupled integro-differential equations, as is 
the case in reference (7). Finally, Sec. IV gives the results of numeri- 
cal integration and their interpretation. 

The numerical integrations were carried out for all partial waves, 
where in higher partial waves Born approximation were used. The 
transition between the eigenstate expansion calculation and Born ap- 
proximation takes place when the results of the two calculations agree 
closely. 
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II. FORMULATION 


A. Derivation of the Differential Equations 
Since spin orbit interaction of the electrons are neglected, the total 
orbital angular momentum L, and the total spin angular momentum S, 
are separately conserved. We can then divide the interactions into anti- 
parallel spin states where S = 0 , and parallel spin states where S = 1. 

We then deal with spatial wave functions of the electrons only and for 
brevity we call the orbital angular momentum the angular momentum. 

Neglecting the motion of the proton of the hydrogen atom and taking 
its position as the origin of the coordinate system, the Schroedinger 
equation for the system can be written 

[H-E] 0(r r r 2 ) = 0, (2.1) 

where r x and r 2 are the position vectors of the bound and free electrons, 
and in atomic units 


H - E = 




- + — -E, 


( 2 . 2 ) 


where E is the total energy of the system andr J2 is the distance between the 


two electrons. We expand the total wave function \p ( r x , r 2 ) in terms of the 


eigenfunctions of the total angular momentum L, 


'M r i> r 2 ) = ^ ^ L ( r i’ r 2 )- (2. 3) 

L=0 
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Since these eigenfunctions are orthogonal and distinct, substitution of 
Eq. (2.3) in Eq. (2. 1) gives. 


[H-E] ^ L (r r r a ) = 0. 


(2.4) 


The explicit form of 'P L (t v r 2 ) is given by 


^( r l’ r a) =( 1 + #\a) H L 

n ”l m 2 1 2 


r 2 ' lu ( k n/ 2 ’ r 2 ) Y ^ 2 o , 2 ( Q 2 ) 


(2. 5) 


<P( n A m l> r i) = r i lp ( n i^l - r i> (2,6 ^ 

Here <p(n 1 'i 1 m 1 ,r 1 ) is the hydrogen atom wave function with radial part 
r' 1 P(n 1 ^ 1 , rj) and angular part Y£ m (fij) and quantum numbers ; 

r- lu ( k n^ 2) r 2 ) is the radial part and m (^ 2 ) is the angular part of the free 
electron wave function with quantum numbers k n ^ 2 m 2 . The relationship 
between the wave number k and is given by 


k 2 = 2 E + 



(2.7) 


l 

Finally the constants C m 1 m 2 M =(^ 1 ^ 2 m 1 m 2 | LM) are vector coupling coefficients 
which make the linear combination of the products of the one electron wave 


functions in Eq. (2. 5) the eigenfunction ofL. In the problem under con- 
sideration n t = 1, 2; \ - 0,1; l 2 = |L -4, | |L+^ | ; mi = . . .,l l 

and m 2 = - / t 2 , . . . , . To make the total wave function symmetric for 
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antiparallel spins or antisymmetric for parallel spins, the operator P 12 
interchanges rj and r 2 while /3 is +1 for the first case and is _1 for the 
second. 

By taking L perpendicular to the z-axis M = 0 and m 2 =-m,. Eq. (2.5) 
can then be written 

'M'l’fj) =(1 +^p 12 ) 21 H c p( n i^i m i> r i) 

n x l x l 2 m l 1 1 

* r 2* u ( k ni V r 2 ) Y -e 2 m 2 ( fi 2 )- ( 2 -8) 


In order that ^P L (r v r 2 ) closely approximates the exact wave function, 
we minimize the expectation value of the energy operator with respect to 
the radial parts of the free electron wave functions, 

8 J'/'l ( r r r 2 ) [H“ E 3 V' L ( r i' r 2 ) d3f i d3r 2 = 0- (2. 9) 

It has been shown by Kohn 10 that the differences between the scattering 
amplitudes obtained from these equations and the exact scattering ampli- 
tudes are quadratic in the difference between t// L (rj, r 2 ) and the exact wave 
function. When the variation is carried out inside the integral we obtain 


J2 C ^^io| ( P*( n l' e i m r r l) y J 2 m 2 ( n 2) [H-E]<// L (r r r 2 ) d 3 r 2 ,dft 2 = 0. (2. 10) 
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By means of Eqs. (2. 2), (2. 8), the Schroedinger equation for the hydro- 


gen atom, 




?(n 1 't 1 m 1 ,r 1 ) = ~ , r x ) , 


( 2 . 11 ) 


and Eq. (2. 7), Eq. (2. 10) reduces to 


E E E C "r"i° C "jV° /?*< n iVr t i> Y V,<n 2 >( l+/8Pia> 

§ ip I / 1 

m 1 n 1 ^ 2 m 1 




+ k’. + 2|f - i 

r 2 r l 2 


<P(n 1 'V m i'- r i) 


x r j 1 u(k , , r 2 ) , , (fi 2 ) d 3 fi dfij - 0 , 


( 2 . 12 ) 


2 2 


where V 2 is the radial part of V 2 By orthogonality of the hydrogen atom 
r 2 2 

and spherical harmonics wave functions, the relation 11 




(2. 13) 


the integration by parts of the exchange terms, and the relation 12 


l 2 'l 1 L _ L-'fci c ^i^ 2 h 

-mjmjO \ / m 1 " m 1 © * 


(2. 14) 
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Eq. (2. 12) leads to 


„ 4 4 + 1 , 2 

V 2 — + k 2 + -2 

r 2 r 2 n x r 2 


u ^2 ’ r 2^ 


-E EE ‘‘A 4% S 

m l n l ^1 4 m l 


’ r l) Y V 2 ( n 2 ) 


l 12 


(fin/ ly ttly, ty) 


U ( k n, ^2 ’ r 2 ^ 


Y 4m'(fi 2 ) +/3T(ni^m;,r 2 ) 


2 2 


u(k ' i’, r,) 




d s r, dfij 


/8 2] + J '(".V 'l> r 2> 




U (k ni< = 0 . (2. 15) 


If 1/r !2 is expanded in terms of the Legendre polynomials and use is made 
of the addition theorem 13 we obtain 


°° \ 

rTT 2Z Pk(COS e,l) 

0 > 


00 +\ \ 

= z] 2Z wk 7XT7 Y x, < n .> 


\=0 


(2. 16) 
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In this expression 0 12 is the angle between the position vectors r x and 
r 2 at the origin, and r< is the smaller and r> is the larger of |r,| and 
l r 2 I . We also introduce 


f r 2 

y^nln'l', r 2 ) = r^O +1 ) P(n4, r x ) P(n ‘ l' , r x ) r^di-j 

-o 

oo 

+ r 2 ^ J P(n^, ri )P(n'r, r 1 )r 1 -<^ 1 >dr 1 . 
r 2 


Then it follows that 


r 


P(ri4, h)P(n'l', r x ) 


00 +X 


12 


‘*-LL 


X.= o u=-X. 


By means of Eq. (2. 18), the relation 14 


i 


Yo Y. Yj, dfi = 

^3 m 3 '“2 m 2 4/ l m l 


'(2Vl)(2Vl)' 


\1K 


4tt(2V1) 


£ ^ 1 ^ 2^3 £ ^ 1 ^ 2^3 
m I m 2 m 3 0 0 0 


and the definition 


(n^|k n 4') = 





P(nt, r) u(k 4', r) Jr , 


Eq. (2. 15) when multiplied by r 2 gives 


(2. 17) 


r 2 ).(2. 18) 


(2. 19) 


( 2 . 20 ) 
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d’ ^VD „ 2 2 

— 7 1 — + k ^. + - 

S* l 1 7 


dr/ 


r 2 ) 


2 (3 

\ U 


•ni L, L.tt «. 

1+1/ mj nij / X=0 /X- - \ 


'2^.' +l\^ XV^l r^'^'fi'2 2^2 t JP * J? ' \ l hr J> / 

-■ 1 M C 1 1 C 1 1 C 2 2 , C . K n.ln, i, ,rJu(k ,^ 2 

2^ # _|_ j 0 0 0 0 0 0 X' 1 * * * * * n l 


. r 2 ) 




(2lJ + lV* 

w^v 


c U ’*- c^‘c U ’V c «A', W v 

AM«2 ®1 0 0 0 /xn^mj' 0 00 x 1 1 


r 2 ) y^l^n 




+/S 


£ 


(-) 


L 


st-t.'-e* • V.> + k 2, .) " (n.' ■ V • r 2> < " A l k „ ,• f 2 ) 


= 0. (2.21) 


In the exchange integrals above we have defined P(k n ^ 2 , r) as u(k n ^ 2 , r) . 
The vectorial equations 

L = i 1 + i 2 , (2.22) 

where L is constant but and i 2 take the values given before, can be 
divided into two groups, one with L ~'t 1 -^ 2 even and the other with 
L - 'tj “^ 2 odd. Since the total spacial wave function has the parity of 
+ ^ 2 , in the first group the wave function has the parity of L and in the 
second a parity opposite to L. By conservation of parity we have two 
distinct groups of interactions. In this problem, where Is, 2s and 2p 
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states of atomic hydrogen are taken into account, it is easy to see that 
when L- is even, the set of quantum numbers k n has 4 

values; one for each of the Is and 2s states and two for the 2p state. 
When L-'fcj - is odd, has one value which corresponds to the 

elastic scattering of electrons by the 2p state of the hydrogen atom. 

Eq. (2.21) is evaluated for these cases and the resulting differential equa- 
tions are listed in Appendix I. In evaluating Eq. (2.21) it should be noted 
'i Z ■i 

that the C 1 2 3 coefficients are subject to the condition that •'C.'t-'C, 
ra l m 2 m 3 1 * 3 

form a closed triangle and m 3 = m x + m 2 . 1S This limits the summation 
over k and /x considerably to few terms only. Summation over 
nij , nij , \ and /x are carried out using the numerical values of the C coef- 
ficients given by Condon and Shortley. 16 

Percival and Seaton 6 have derived the same differential equations 
for the scattering of free electrons by atomic s, p and d electrons in 
hydrogen atom using the theory of irreducible tensor operators to 
evaluate the interaction terms between the two electrons in the dif- 
ferential equations. The calculation becomes considerably simpler in 
this way. The results of the two methods are identical. 

In the rest of the paper except Sec. IIE we discuss the solutions 
to the four coupled differential equations given in Appendix I and which 
arise when L-'t 1 ~^ 2 * s even * The single differential equations for 
L - / C 1 - ^ 2 odd is derived in Sec. IIE. Its numerical integration can be 
treated as a special case of the four coupled differential equations. 
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When the integrals representing the direct potentials in the four 
differential equations are evaluated and some change is made in the 
limits of the exchange potential integrals, these equations can be 
written in the following matrix form 


U ( k n / n* r ) = 2 V U ( k n / n* O' ( 2. 23) 

The four components of u are the four radial functions of the free electron. 
V is a 4x4 symmetric matrix that is the sum of three matrices. 


d> , , 


dr ‘ 


+ k - • 


V.. = D. . + E.., 
i j i j i j 


v 1 r°° 

E. . = F. . + / g v . h v . dr, 

v-\ •'O 


(2.24) 


where is the direct and E £ . is the exchange potentials and both are 
functions of r. The matrix E.. contains in addition integrals with respect 
to r and for the purpose of numerical integration it can be written as the 
sum of two matrices. The explicit forms of D^, F i *, g and hV. are given 
in Appendix II. The value of a is Z for i=j=3 and i=j=4, and is 1 for all 
other values of i and j. It is understood that for the exchange terms 
the components of u on the right hand side of Eq. (2. 23) are inside the 
integrals of the exchange terms. 
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B. Derivation of the Transmission Matrix From Solutions 
of the Differential Equations 

The method is similar to that used by Bransden and McKee, 17 
and by Marriot 4 Eq. (2. 23) constitutes a set of 4 coupled, second 
order, differential equations. Three components of u can be elimi- 
nated from these equations, resulting in an 8 tn order differential 
equation for the remaining component. Therefore there are 8 sets of 
solutions to Eq. (2.23), However, only half of these solutions are 
regular at the origin. Each of the four regular solutions corresponds 
to a definite vector u. The four vectors can properly be represented 
by a 4x4 matrix u nj ,n,j = 1, 2, 3, 4, where n corresponds to the partic- 
ular component and j corresponds to the particular solution of u. The 
four solutions are carried out numerically in the next section. 

From the explicit form of V it can be seen that V vanishes at infinity. 
The asymptotic solutions of u as given by Eq. (2. 23) is therefore 

U nj ( r ) ^ a nj sin ( k n r - + 8 nj )> (2.25) 

n, j - 1, 2, 3, 4, 
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where a nJ is the amplitude and S nj is the phase shift of the j th solution of 
the n th component of u. 

Corresponding to the 4 components of u there are 4 channels open to 
the reaction. If the incident wave be in the m th channel (m =1, 2, 3, 4), the 
travelling wave in the n* channel will be given by 


u n( r ) ~ exp 




8(m, n) 


s m „ ex P 


H(k„r 



n = 1, 2, 3, 4. 


(2. 26) 


The constants S mn are the amplitudes of the scattered waves. Since Eqs. 
(2.26) are also the asymptotic solutions of Eqs. (2.23), they must be 
equal to linear combinations of Eqs. (2. 25). If we call the coefficients 
of the linear combinations we must have 


4 



n, m = 1, 2, 3, 4. 


(2. 27) 


4 


On the right-hand side, we have used the normalization of Blatt and 
Weisskopf. 18 If we equate the coefficients of exp£-i^k n r -i / n 77^j ancl 
exp [‘( k " r ~5 ( "’ 7 )] in Eqs. (2. 27) we obtain 


7~! P i 3 nj eX P “p S nj] = 8 (m.n) 


j=l 


L 

j=i 


-2i 


P. a nj exp [i S nJ ] = S„ 


Separation of Eqs. (2.28) into real and imaginary parts gives 


4 

[( Rp j) sin S nj -(A Pj) cos S nj ] a nj = ~ S(m,n; 
l-i * " 

4 

Z [ (BPi) cos S n . -h PjJsin a n . = 0 


j 


4 

[( RP i) Sin 8 n, +(^ P j) COS S nj] a „j = 5= 

J.l Y " 


R S_ 


2^ [(RP,) COS S nJ -(A Pj) sin § nj ] a nj = -~ A S, 
j =1 f n 


(2. 28) 


(2.29) 


In the above R or A represent the real or the imaginary part of the quantity 
that follows them. Eqs. (2. 29) are a set of 16 linear equations for 16 
unknowns RP. , A P., R S and i S Once these unknowns are found, 19 the 

j 9 j 7 mn ntn 7 

magnitude of S mn will be given by 
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|s mn | 2 = (RS mn ) J + (c9 S ran ) 2 . 


(2. 30) 


The cross section is obtained by asymptotic expansion in spherical 
harmonics of the incident plane wave 20 


exp [ikz] 


(2/ + 1) M i /+1 
l = 0 




“ exp 



t , o 


(2.31) 


The magnitude of the ingoing wave on the right hand side of Eqs. (2. 27) for n = m 
is [k m /77(2/ m + 1)] 1/2 times the magnitude of the partial wave of the ex- 
pansion ofr exp [ i k m z ] . . The plane wave has a flux of v which, in atomic 
units, is equal to k . The ingoing flux of the right-hand side of Eqs. (2. 27) 
is therefore k 2 / [rr( 21 +1}]. The outgoing flux in the channel n i- m is 

m ' ra ' ' 

|S mn | 2 . The cross section is obtained when we average the ratio of the 
outgoing flux to the ingoing flux over the initial states, and sum over the 
final states. For a particular spin state of the two electrons, unpolar- 
ized electron beam and unoriented atoms, the multiplicity of the initial 
states is (2/j + 1)(2 l 2 + 1) , where and l 2 are the angular momentum 
of the bound and free electrons. For a polarized beam, m 2 = 0 , where 
m 2 is the magnetic quantum number of the free electron. Then m 1 - M » 
where m 1 and M are the bound electron and the total magnetic quantum 
numbers. Since M is constant there is only one initial state for a 
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polarized beam. The multiplicity of the final states is 2L+ 1, where L is 
the total orbital angular momentum. Since the cross section for 

m 2 

m /n is 


Q 

^ mn 


7T ( 2L + 1) | S 

(2/j 


ra ^ n . 


(2. 32) 


The outgoing partial wave in the incident channel m consists of the 
scattered wave plus the outgoing wave given in the expansion of the plane 
wave. Then, according to Eqs. (2.27) for n = m the magnitude of the 
amplitude of the scattered wave is |l-S mm |. The elastic scattering 
cross section is therefore given by 


Q 

mm 


7T (2L+ l) 

k^(2/ 1+ l) 


1 - S. 


(2.33) 


If we define a matrix T by the relation 

T = 1 - S , (2.34) 

Eqs. (2. 32) and (2. 33) can then be combined into a single equation. 


Q..= ”< 2t - +1 > I T J 1 - (2.35) 

ki(2i,+ l) 

X is the transmitted amplitude in the n^ channel due to an incident wave 

mn 

in the III th channel. The elements of T „ constitute the transmission 

mn 

matrix. 
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The matrix S has two properties that are useful as tests on the 


accuracy of numerical integration. From Eq. (2. 26) it can be seen that 
S transforms the ingoing wave into the outgoing waves. The continuity 
of the electronic current requires that S be a unitary matrix 

I S mn | 2 = 1, m= 1, 2, 3. 4. (2.36) 

nsl 

Furthermore, since the Hamiltonian is Hermetian, S must be sym- 
metric, 18 


S 


ran 


= S. 


(2.37) 


Eqs. (2. 36) and (2. 37) are used as tests on the accuracy of numerical 
integration. 


C. A Useful Relationship 

A relationship based on the symmetry of the interaction potentials 
can be derived which serves as another test on the accuracy of the 
solutions. The / th and the k th solutions of the i th component of u by 
Eq. (2. 23) are given by 


i 
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+ kj - li 


0 






d J 




/ «i + *>“ 

i ;; 


u ik 



u 


jk 


j 


(2.38) 


Multiplying the first by u ik and the second by u kl , subtracting the two 
expressions and summing over i gives 



(2. 39) 


Since V ij = V jl » the interchange of the summation indices changes the 
sign on the right hand side of the equation, the right hand side must there- 
fore be zero. Integrating the left hand side from zero to infinity we obtain 



= 0 . 


(2. 40) 


Integrating the above equation by parts, and applying Eq. (2. 25) we obtain 


4 

L 

lal 


l ik 3 i / 


a il Sil1 ( 8 ik “ 8 i/) = 0 


k, / = 1, 2, 3, 4, Vf l . 


(2.41) 
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Although the terms containing the exchange potentials do not cancel out 
on the right hand side of Eq. (2. 39) the cancellation does take place after 
the integration is carried out in Eq. (2.40). 


D. Transmission Matrix According 
to Born Approximation 

The Born approximation consists of neglecting the exchange potential 
terms appearing in the V matrix of Eq. (2. 23), and also of neglecting 
all the direct potential terms in this matrix except those terms that 
connect the incident channel to all other channels. 21 Eq. (2. 23), when 
the incident wave is in the m^' channel, reduces to 



('n + 1)' 


“n = 2 D nm u, 


(2. 42) 


n = 1, 2, 3, 4, 


u m and u n are given asymptotically by 


u m ~ sin ( k m r ~ l m 7T /' 2 ) (2.43) 

u n ~ B nm cos (k„r - l n n/ 2). (2.44) 

We have chosen the constants of proportionality of and u n such that 
B nm is the Born approximation of the reactance matrixR. 22,23 Eq. (2.43) 
shows that u m must have the following form 24 
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U m - k m r J/ m ( k ra r )‘ 


(2.45) 


where i, (k r') are spherical Bessel functions. Furthermore, if y repre- 
U m K m ' n 

sents the homogenous solution of Eq .(2. 42), it must have the following 
forms 

y n = a n k n r j/„ ( k n r ). ( 2 * 46 ) 

y n - 3 n Sin ( k n r - l n rr/2), (2.47) 

with a n some unknown constant. Multiplying Eq. (Z. 42) on the left by 
y n and integrating the result from zero to infinity we obtain by partial 
integration 


r 

1 y D u dr = 

r 00 

y n 

rtf 

n 5 

t J 

> 

Ldr 2 

r J 


r~ 


00 


r a 

d 1 


— 

V — u - u — y 
* n dr n n dr n 





0 

— 

-k M 

n 

a E . 

n n m 



u dr 


The last equality has been obtained by noticing that y n and u n vanish at the 
origin, and by using their asymptotic forms as given by Eqs. (2.44) and 
(2.47). We therefore have 


B = - 2(k k ) 

nm ^ n m / 


| n ( k n r ) D nn 


Urn ( k m r > ^ dr - 


(2.48) 


This is identical to the expression given for B by Seaton. 
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The transmission and the reactance matrices are related by T = - 2i R/(l - iR ) . 
Since B « 1 , the transmission matrix according to the Born approximation is 
given by 

T n m B = 4i ( k n k m > 1/2 j Un ( k n D nm Urn ( k m r > r2dr ' ( 2 - 49 ) 

•'o 

Substitution of Eq. (2.49) is Eq. (2.35) would give the cross section according 

Si *' 

to the Born approximation. It should be noted that the symmetry of T insures 
Eq. (2.37) to be satisfied while Eq. (2.36) is not satisfied anymore. If we 
define T = -2 iB/(l - i B ) , then S = (1 + iB )/(l - i B) and the symmetry of B 
makes S unitary. Eq. (2.36), the continuity of current, is then satisfied. 

This is the definition of approximation II of the Born approximation according 
to Seaton 23 which will be discussed later. 


E. Elastic Scattering of Electrons 
by the 2p States of the Hydrogen Atom 

The angular momentum of the free electron, l 2 , in the 2p channel has 

the values L- 1, L ,L + 1 , where L is the total angular momentum of the 

system. The first and the last values were considered in previous sections. 

The case l 2 = L corresponds to a wave function in the 2p channel with a 

parity different from all channel wave functions considered previously. 

It therefore corresponds to elastic scattering. The wave function in this 

case is given by 
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(2. 50) 


Tl 


= <» * 0P„>2] u(k ‘ t L,,i> V.,< n a >- 


When Eq. (2.9) is formed with this wave function, and minimized with 
respect to u(k 2 L,r 2 ), treatment which led to the derivation of the four 
differential equations will give the following differential equation 



_ L(L + 1) + 2 
r 2 r. 


u L ( r ) 


+ S(L,1) r R 21 (r) (2p|k 2 L) 

+ 2 |y 0 (2p2p, r) --|-y 2 (2p2p, r)J u L (r) 

= 0 (2.51) 


+ 2ySrR 2I (r) 


3 y L -i ( 2 P k 2 L > r ) 3 y L +i ( 2 P k 2 L > r ) 
r (2L+1) (2L-1) + (2L+1) ( 2 L + 3) J 


The asymptotic solution of this equation is given by 


u L ~a L sin(k 2 r- l/n/2 + S L ) . (2. 52) 

If the scattering amplitude is designated by T 5S it can be shown from 
Sec. IIB that for a particular L 

T 55 = 1 - exp 2i 8 = - 2 i exp (i8)sin8, (2. 53) 

where for simplicity we have suppressed the subscript!,. The cross 

section, according to Eq. (2.35), is given by 

O = 4tt( 2L+ 1) sin 2g (2.54) 

3 k 2 
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The total elastic scattering cross section by the 2p states is the sum of 
this cross section and the cross section corresponding to l 2 = L - 1 and 
l 2 = L+ 1 which were considered previously. 

The Born amplitude, Eq. (2. 48) in this case is given by 

B 55=-2k 2 f j L (k 2 r)D ss j L (k 2 r)r 2 dr, (2.55) 

0 

where by Eq. (2* 51), 

D ss = “ 7 + y 0 (2p2p> r ) - y 2 (2p2p. r ) • (2.56) 

III. NUMERICAL INTEGRATION 

A. Decomposition of the Differential Equations 
If it were not for the definite integrals appearing in the potential 
matrixV, the set of the four coupled differential equations (2.23) could be 
integrated by any standard technique. The presence of these unknown 
constants whose integrand involve the unknown functions makes it 
necessary to solve these equations by iteration or by transformation 
of u into other vectors, whose differential equations do not contain 
definite integrals. Since the terms containing definite integrals are 
small as compared to the direct potentials, the iteration method can be 
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used by assuming that the values of these integrals are zero. The 
differential equations are then integrated, the values of the definite 
integrals that are subsequently obtained are substituted in the differential 
equations, and the integration is repeated. The process is repeated 
until sufficiently consistent values of these integrals are obtained. This 
method is useful if the convergences of the constants are fast enough, 
and the cross section is not very sensitive to the values of these constants. 

In the second method, the transformation of u fixes the values of the 
constants and thus avoids iteration, whereby the computation is reduced 
considerably. The description of the method will be given here. 26,27 

By making use of Eq. (2. 24), Eq. (2. 23) can be written 



M'i+iy 


u. = 


2 



+ F. 


ii> u i 



(3. 1) 


where 



(r) dr. 


(3.2) 


We introduce the functions v £ and that are solutions of the following 
differential equations 
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"d 2 

+ k 2 - 

h 0 i + 1)' 

< 

11 

K3 

dr 2 

1 

r 2 

l 

"d 2 

+ k 2 - 

l i <*i + 1)1 

1 1 k/ — 

dr 2 

r 2 

u i - 


. = 2 £i 

j=i 


u ;' = 2 21! i?‘i +F ‘i]< 

j »* 


(3.3) 


+ 2 S (i, k) gjj. 

Then u. is given by the following expression 


(3.4) 


u. = v. + 


k = 1 /si Vs 1 


kl 


(3.5) 


Eq. (3. 5) can be verified by multiplying Eq. (3. 4) by C£/ , summing over 
k, / and v, and adding to Eq. (3.3), where upon Eq. (3. 1) results. 
Substitution of Eq. (3. 5) in Eq. (3. 2) gives 


4 4 <7 


[S(ij/*. Wvj-Aff'] q, =Bf, , (3.6) 

k= 1 /si V«1 

i, j =1,2, 3, 4; /x = 1, 2 for i = j = 3 and i = j = 4; ^ = 1 otherwise, 
where and B^j are defined by 


a?:' 


= { h fi 


dr 


*-f 


= | h* v. dr 


> * 


(3. 7) 
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The numerical integration is carried out by integrating Eqs. (3. 3) 
and (3. 4) by any standard method, calculating and B^. by Eqs. (3. 7) 

and, finally, solving the system of 18 algebraic equations given by Eqs. 

(3. 6) to find . With the known values of these constants the integra- 
tion of Eqs. ( 3.1 ) is straightforward. 

The determinant of Eqs. (3. 6) becomes singular for L = 0 and 1. This 
is shown in Appendix III. To remove the singularity, some of the Cjj are 
chosen arbitrarily, and the rest of the Cjj are found in terms of the chosen 
ones. 


B. Solution at the Origin 

In order that the four soltuions of u be independent of each other we 
must have 


7^ C. u,. /0, i = 1, 2, 3, 4, (3. 8) 

jal 

where C ^ are some constants. A necessary condition for this to be 
satisified is that the determinant of Eq. (3. 8) be nonzero. 


1 1 1 1 f 0 - 


(3.9) 


It is not difficult to see that this also is a sufficient condition. At the 
origin the solution u^. can be expressed as power series in r, 


00 



v-0 


s. 


“u 


(3. 10) 
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where a^. are the coefficients of expansion, and are given integers 
for each component of u and are fixed by the behavior of Eq. (2. 23) 
at the origin. We can satisfy Eq. (3. 9) at the origin by having 


1 1 a Ij 1 1 / 0. (3.11) 

By choosing suitable values of a°., subject to the restriction (3. 11), 
four independent solutions are obtained. 

C. Solution at Large r 

With given initial values the solution of Eq. {2. 23) can be extended 
from origin to any desired value of r. In order to obtain the 
asymptotic amplitudes and the phase shifts, the presence of the 
centrifugal and long range potentials make it necessary to extend 
the solutions to infinity. This is undesirable because of the time 
consumption on the computer, and the accumulated error due to 
the long range integration. Seaton 23 has solved the problem of 
r“ 2 long range potentials occuring in the off diagonal terms of the 
potential matrix V by diagonalizing the asymptotic form of the dif- 
ferential equations (2. 23) and the corresponding S matrix. By an 
inverse transformation the elements of the original S -matrix are 
found. 

Instead, we develop here a perturbation theory which is based 
on the method described by Mott and Massey. 28 The error in the 
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resulting solution is inversely proportional to the squared of the 


distance from the origin. 

Eq. (2. 23) for large distances of r can be written 


r ^2 


Ldr : 


+ k. 


«(k n / n , r) = 2 Uu(k n r), 


(3. 12) 


where U is the sum of the centrifugal potential matrix and the asymptotic 
form of the V matrix. The elements of U are given in Appendix IV. 

A component of Eq. (3. 12) is of the following form: 


-^-+k 2 ] u(r) = g(r) 

Ldr 2 J 


g(r) « k 2 u(r) , g(r ) -0 as r -°o 


(3. 13) 


The perturbation theory is applied between some large distance R and 
infinity. Suppose u vanishes atR, then we have the following boundary 
condition 

u(R) = 0 . (3.14) 


If we represent the solution of the homogenous equation by y(r), at 
infinity we must have 


y(r) = a s in (k r - kR) I 

u (r) = ( a + Aa ) sin (k r - k R + 77 ) 


(3. 15) 


where AA and 77 are generated by g(r) . Since g(r) is small, we can write 
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u = y(l + O . 


(3. 16) 


where 4 is a small function. Substitution of Eq. (3. 16) in Eq. (3. 13) gives 



y, 


where, upon double integration, we obtain 


£ = 



g(r') y dr'. 


(3. 17) 


(3. 18) 


The constants of integrations are fixed by the condition (3. 14) and the 
fact that u' (R) = y'(R> 

We now integrate Eq. (3. 18) by parts, 


£ = 



g(r) y 




g(r)ydr 

‘'r 



(3. 19) 


When the integration with respect to y is carried out, and the result is 
substituted in Eq. (3.17), we obtain 


u(r) = sin (kr - kR) 


hr 


g(r) cos (kr - kR) dr 


(3. 20) 


[-if 


+ cos (kr - kR) --j- I g(r) s in (kr - kR) dr 
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Comparison of the second of Eqs. (3. 15) and Eq. (3. 20) shows that 


■if 




Aa = -f- | g(r) cos (kr - kR) dr 


> » 




(3. 21) 


g(r) sin (kr - kR) dr 
r J 


to first order. The functions g(r) in the four differential equations (3,12) 


are given by 


g l (r) = 2 ^U li u ) . 


(3. 22) 


To first order this can be written by 


gi(r) = 2 


L 


u 


ij 


sin (k. r - k, R.) 


(3.23) 


where Rj is the last zero of with positive slope. Substitution of this 
equation in Eq. (3. 21) gives 


A a 


r Ri 

.-Zx I 




cos (kj r - k. R^U. . sin(kjr-k.R. )dr 


y. (3.24) 


r Ri 

7, = 23 A. J sin < k .*- - k i R i> u ii sindtjr-k.Rpd, 


A a and 77 . can easily be calculated by substituting the values of 
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u. . 

1 J 


from Appendix IV, integrating the resulting integrals by parts and 


retaining the leading terms. 

The asymptotic amplitudes and phase shifts are given by 

a. (®) = a. (R.) + Aa. "1 


> 


(3.25) 


S.(oo) = 3. (R.) + 7). + [L- 8(i,3) + 8(i,4)]7r/2 

J 

where a.(R.) and S.(R.) are the amplitudes and total phase shifts calcu- 
lated at R. by the machine, and where S(i,3)and S(i,4) are the S 
functions. 


D. Details of the Numerical Integration 
Milne 1 s 29 method with variable mesh size and Simpson's 30 rule were 
used for the integration of the differential equations and evaluations of the 
integrals respectively. As the solution advances from the origin, the 
differential equations become less sensitive to the size of the increment, 
and the error of integration falls below certain small number e. At each 
value of r the value of the function is found, first with the given value of 
the increment, and second with the value of increment divided in half. 

The error of integration is defined as the difference between these two 
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solutions. When the error becomes small the increment is doubled until 


a maximum value is reached. At some distance Rj all the exchange 
potentials and, similarly, all the direct potentials except those represent- 
ing optically allowed transitions and the 2p-2p elastic scattering potential 
become vanishingly small. (See Appendix IV. ) At this distance the 
set of differential equations is replaced by the simpler set con- 
taining only these potentials. The integration is continued until some 
distance R 2 , where the first order solution of the rest of the range of 
integration is obtained by the method developed in Sec. IIIC. No 
attempt was made to solve any set of linear equations or any matrix 
equations, as these equations are solvable by the computer in their 
original form. 

The values of the constants of the numerical integration are given 
below. h. and h f are the initial and the final increment of integration. 

In some exceptional cases, different values were used. 

Table III. 1 
h. 1 x l(T s 

h f 0.05 

e lx io -4 

Rj 30 

R 2 200 

All quantities are in units of Bohr radius except e which is dimensionless. 
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rv. RESULTS 


The four differential equations listed in Appendix I 
were integrated by the methods described in Sec. III. Based on the 
Eqs. (2.41), (2.37) and (2.36) we define the three relationships, 

4 

/ k. a. a s n sin(S. - 8. ) 

/ i im in v lm in' 

D =J-1 , 

mn 4 

y k. a. a. |sin(8. -8. )| 

/ 1 im in | v 1m in'| 

i » 1 

m,n = 1, 2, 3, 4, m i n, (4.1) 


D'_ 



l S „ 


+ S 


m, n = 1, 2, 3, 4, m ^ n, 


(4.2) 


D 


n 

m 



n c 1 


m = 1, 2, 3, 4. 


(4.3) 


In an exact solution of the four differential equations the right-hand side 
of these equations should vanish, they can therefore be used to test the 
accuracy of the numerical integration. As an illustration the numerical 
values of D mn , D' mn and for the case of ls-2s-2p eigenstates coupling, /3 = +1, 
k x = 2.0 and L = 3 will be given below 
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D ia 

= 1.4 x 10- 3 , 

D 13 

= 2.6 x 10- 4 , 

D 14 = 1.3 x IO' 3 

» 23 

= 5.1 x 10* 4 , 

D 24 

= 2.2 x 10* 3 , 

D 34 = 1.8 x 10* 3 

D 12 

= 7.6 x 10- 4 , 

*13 

= 5.1 x 10* 3 , 

D' 14 = 5.6 x IO' 3 

d' 23 

= 5.4 x 10* 3 , 

D 24 

= 5.7 x io* 3 , 

D' 34 = 1.3 x 10* 3 


D* = 1.8 x 10- 4 , D" = 7.6 x 10* S , D" = 2. 5 x IO' 4 , D" = 4.4 x io* 6 . 

When L is large, particularly at the threshold, the cross sections are 
small and the values of these cross sections are more sensitive to the 
errors in the numerical integration. As a result the D values become 
large. In the tables that follow, except for a few cases, the cross 
sections that are listed have maximum D values less than 1%. In the 
exceptional cases by combining different independent sets of solutions 
and different mesh sizes it has been assured that the cross sections 
listed are accurate to within a few units in their last significant 
figures. 

The cross sections listed in the tables are obtained by averaging 
the values of cross sections from two independent sets of solutions of 
the four differential equations. These independent solutions are 
obtained by choosing two different values for the determinant (3.11). 
For small values of k x a maximum mesh size of 0.1 in Bohr radius is 
used while for large values of kj this mesh size is 0.05. It is believed 
that within the framework of the present formalism the cross sections 
are accurate within one or two units in their last significant figures. 
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Corresponding to the three channels Is, 2s and 2p, the cross 
section matrix is a 3 x 3 matrix. By Eqs. (2.32) and (2.37) the excita- 
tion cross sections 2s — »ls, 2p— -»ls and 2p — »2s can easily be found in 
terms of the excitation cross sections Is — >2s, Is— >2p and 2s — »2p. 

The excitation cross section 2s — >2p is quite large and should be calcu- 
lated by taking into account the energy difference between the levels 
2Sj /2 , 2P 1/2 and 2P 3/2 . The summation over partial cross sections does 
not converge for this transition and the method of partial wave is not 
useful. The Born calculation of this cross section is given by Seaton. 31 
In the 5 tables that follow the excitation cross sections Is - 2s, Is - 2p, 
and the elastic scattering cross sections Is- Is, 2s-2s and 2p-2p are 
given. The ls-2s, ls-2p, andls-2s-2p eigenstates coupling approxima- 
tions correspond respectively to the inclusion of the Is, 2s; ls,2p and 
Is, 2s, 2p hydrogen eigenstates in the total wave function. The case /3 = 0 
corresponds to the neglect of the exchange terms in the ls-2s-2p eigen- 
states coupling. 

To find the contribution of the higher partial waves whose cross 
sections are difficult to calculate by the method of eigenstate expansion, 
we use the regular partial wave Born approximation. Although the 
partial wave Born approximation as defined by Seaton, 23 in which the 
continuity of current is satisfied, agrees better with eigenstate expan- 
sion calculation for particular elements of the cross section matrix, this 
is not true in general. As a result, for the sake of consistency, we calculated 
the contribution of higher partial waves through regular Born approximation. 

The cross sections corresponding to the 5 tables are plotted in the 
5 figures. 
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V. DISCUSSION 


(i) The results of numerical integration using noniterative method 
agrees with the iterative numerical integration of Burke, Smith and 
Schey, 7,32 although the two results are different sometimes by as much 
as 10%. The difference between the two methods is large when the 
cross section is small. 

(ii) ls-2s excitation cross section. As is seen from Fig. 1, the 

Is -2s and the ls-2s-2p eigenstates calculations are in good agreement 
with the experimental results of Lichten and Schultz 1 although the 
calculations do not give the sharp peak in the cross section at the 
threshold as is observed experimentally. The shape of the exchange 
neglected ls-2s-2p eigenstates calculation curve is in excellent agree- 
ment with the shape of the experimental curve. The two curves coin- 
cide beyond 25 ev. We notice that the Born approximation agrees with 
the eigenstate expansion calculations for values of bombarding 
energies close to the threshold. The normalization of the experimental 
results in reference (1) to the Born approximation at 45 ev appears 
therefore to be justified. 

(iii) Is -2p excitation cross section. The experimental curve 3 
agrees in shape with the calculated curves but it is lower than all of 
them (cf. Fig. 2). It should be noted that the inclusion of the 2s state 
in the Is - 2p excitation cross section calculation increases this 
cross section. 

Since the experimental curve of Stebbings et al. in Fig. 1 is 
normalized according to the experimental curve in Fig. 2, if in a new 
measurement of the ls-2p transition cross section higher values for 
the cross section is found, the two experimental curves in Fig. 1 will 
be brought into better agreement. 
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(iv) 2s -2s elastic cross section. The cross section at the thresh- 
old is large. A Ramsauer effect is evident in the L = 0 case (cf. Table 
IV). 

(v) 2p-2p elastic cross section. The cross section at the threshold 
is considerably larger than 2s -2s elastic cross section. The Born 
approximation agrees with the eigenstate expansion calculations at higher 
energies. This cross section is not calculated in reference 32. 

(vi) Although we have neglected the spin orbit interaction we have 
investigated the role that spin plays in electron atom collision. Assum- 
ing that electrons are identical but spinless particles, the total wave 
function should be symmetric with respect to the coordinate interchange 
of the two electrons. The curves Is - 2s -2p symmetric in the 5 figures 
correspond to such a wave function. The disagreement between theory 
and experiment is an indication of the important role the antisym- 
metry property of the total electron wave function plays in the free 
bound electron collisions. 

(vii) A study of the 5 figures shows that the Born approximation, 
despite its simplicity, agrees with eigenstate expansion calculations for 
low and close to the thresholds bombarding electron energies. Of 
particular interest is the fact, yet unexplained, that the Born approxi- 
mation underestimates some cross sections and overestimates others. 

VI. CONCLUSION 

The hydrogen eigenstate expansion of the total wave function of 
the system of a free electron in collision with a hydrogen atom seems 
to be the most natural expansion in any partial wave calculation. The 
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formalism can easily be extended to the electron impact atomic hydro- 
gen excitation cross section calculation higher than n = 2, and to the 
electron impact excitation cross section of any hydrogenic atom 
whose central potential field is given by a potential such as Hartree or 
Hartree Fock potentials. Where the Born approximation does not give 
accurate results the eigenstate expansion seems at present to be the 
most suitable method. 

In addition the treatment of the angular momentum as described 
here is the same for any system of one free and two bound particles 
with radial interacting forces, and the formalism developed here can 
be applied to such problems. 

We conclude that the Born approximation can be used with more 
reliability in electron impact atomic excitation cross sections. 

Finally, it is believed that the noniterative technique employed 
here to solve the radial differential equations will find more 
applications in problems containing exchange integrals. 
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APPENDIX I 


The Four Differential Equations 


= 2y 0 (lsls, r) u(kjL, r) + 2y 0 (ls 2s, r)u(k 2 L, r) 


[a(2L + 1)] 


( Is 2p, r) u(k 2 L- 1 , r) 


[s(2L + 1 >] 


(Is 2p, r) u(k 2 L + 1, r) 


2 B O ft 

+ 2171 rR io( r )y L ( lsk i L, r) + 2 LTI rR 20 (r) y L ( Is k 2 L, r) 


(2L+ 1) (2L-1) 


rR 2 i( r )y L1 (lsk 2 L-l, r) 


2/3 3 ( L + *) 

_(2L+ 1) (2L + 


1 ) 

+ 3 ) 2 _ 


rR 2i( r )y L+ 1 ( lsk 2 L +l. r) 


-/3S(L, 0) ( 1 + k, 2 ) (Is |kj L) rR 10 (r) 


-/3S(L, 0)(l + k 2 2 )(ls|k 2 L)rR 20 (r) 


-/3S(L-l,0)(l+k 2 2 )(ls|k 2 L-l)rR 21 (r). 
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±\ki + i 


dr 2 


2 r 


u (k 2 L, r) 


= 2y 0 (Is 2s, r) u(k , L, r) + 2y 0 (2s2s, r)u(k 2 L, r) 


I 1/2 


+ 2 


|3( 2L +1) 


y, (2s 2p, r) u(k 2 L - 1 , r) 


- 2 


T 1 1 1/2 

[ _3 (~ 2L + 1 ) J yi( 2s2 P< r)u(k a L+l, r) 


+ 2L + 1 rR io( r )y L (2sk i L > r) + 2LTi rR 2° (r) y L ( 2sk 2 L * r > 


+ 2/3 


3L 


(2L + 1) (2L - l) 2 


1/2 


rR 


2 i ( r ) y L _ i ( 2s k 2 L - 1 , r) 


- 2/3 


3 ( L + 1 ) 


1 1/2 


(2L + 1) (2L + 3) : 


rR 2i( r >y , +1 (2s k 2 L + l, r) 


-/3S(L,0) (± + k, 2 ) (2s|kjL)rR J0 (r) 


-yS8(L,0) (| + k 2 2 ) (2s |k 2 L)rR 20 (r) 


/38 (L - 1 , 0) (^ + k 2 2 ) (2s|k 2 L-l)rR 21 (r) 


(A2) 
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(L-l)L 2~ 
r 2 r . 


u (k 2 L- 1, r) 


= 2 


1/2 


3( 2 L + 1 ) Y: (ls2p,r)u(k 1 L,r) + 2 3(2L + 1) 


3 


1/2 


+ 2 


L- 1 


y 0 (2p2p, r) + g ^ 2L+1 ^ y 2 (2p2p, r) 


u (k 2 L - 1 , r) 


-§ y2 ^ 2p2p ’ r > u (k 2 L +i, r) 


+ 2/3 


3L 


1/2 


(2L + 1)(2L-1) 2 


rR io( r ) y L .i (2pk, L, r) 


+ 2/3 


3L 


(2L + 1) (2L-1) 2 


1/2 


rR 2 o( r ) y L . 1 (2pk 2 L, r) 


+ 


6/3 

2L- 1 


rR 2 i(r) 


y L (2pk 2 L-l, r) 
( 2 L + 1 ) 2 


L-l 
+ 2L- 3 


y L-2 ( 2 P k 2 L “ 1 - r ) 


6/8VL(L+l) (2pk L + l, r) r R 21 (r) 
(2L+ 1 ) 2 L 


- (L, 1) + kj 2 ^ (2p | kjL) rR 10 (r) 


- /3S (L, 1) + k 2 2 ) ( 2p | k 2 L) rR 20 (r) 

-J3S(L,2) ({ + k 2 2 ) (2p | k 2 L-l) rR 21 (r) . 


i ( 2s 2p , r) u(k 2 L, r) 


(A3) 
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~d*_ 

dr 2 


+ k. 


(L -f 1) (L + 2) 2 

r 2 r . 


u(k 2 L+l 


r) 


= - 2 


- ~ j 1/2 — — 

3(2L + 1) yt(ls2P, r) u(kjL, r) - 2 3 ( 2 L + 1) 


1/2 


VlCL + 1) y 2 (2p2p, r) u(k 2 L-l, r) 
2L+1 


+ 2 


(2p2p, r) + 


(L + 2) y 2 (2p2p, r) 
5( 2L + 1) 


u(k 2 L + 1 , r) 


- 2/3 


3(L + 1) 


(2L+ 1) ( 2L + 3) 2 


1/2 


rR i 0 (r)y L+1 (2pk 1 L, r) 


- 2/3 


3(L-h 1) 


( 2L + 1) ( 2L + 3) 2 


1/2 


rR 2o( r )y L+1 ( 2 P k 2 L - r ) 


OVW'.O 


+ 2l/+ 3 rRzi ^ 


y L ( 2p k 2 L + 1 , r) (L + 2) y L+2 (2pk 2 L + l, r) 

- + 


(2L+iy 


(2L + 5) 


-yS5(L,0) (j+k 2 2 ] (2p| k 2 L + l)rR 21 (r) . 


j (2 s 2P, r) u(k 2 L, r) 


(A4) 
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APPENDIX II 


Elements of the Potential Matrix 
Elements of D^j: 



44 


~1 

_r 


3 r 

+ 4 + 4 + 




, 6(L - 1) 
2L + 1 

, 6(L + 2) 
+ 2L + 1 




1 r r 2 

+ 6 + 24 + 144 
1 r r 2 

+ 6 + 24 + 144 


9 27r 

_+ 64 


y] 


^4=D 41 


/ L + l \ 1/2 [j_ 
\2L + 1/ |y- 


243 


/ 1 3 9 27r\ "rl 

(7 + 27 + 8 + 64j‘ 


D 23 “ D 32 = " 3 ( 2L L + l) [75 " ^ + T + 1 + Z + 54) e ~ f ] 
D 24 = D 42 = 3 (lITi) “(^ + 7 + I + § + 1^) 


D 34 =D 43 =-18 


'L(L+1)V /2 fl (l 1 x _l. 1 _r_ 

_(2l 7I?J P"lr“ + 7 2r + 6 + 24 


Elements of F. 
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833 


3/^1 r L4t f = 
(2L- 1) (2L + 1) 2 ’ “ r L 


g33 


3/3(L-l)R 21 r L - 1 
(2L- 1) (2L-3) ’ h33 



- + k 2 

!L ^— S(L, 2)r 


3/3R 21 r L+1 
(2L + 3) (2L+ l) 2 




- + k 2 
4 + 2 „ 

—~2 S ( L , 0) r 


, 3/3(L + 2)R 21 r L+3 R 21 

g 44 - (2L + 3)(2L + 5) «« ” r L+2 


a / n L+l l n 

812 ~ 2 L + 1 K 2° ’ n i2 R io 


1 l+k 2 2 


-I"— S < L ’°> r 


Lr 


&21 ~ ^12 [*20 ^ * 1 oj > ^*21 “ k 12 [* 10 ~^* 2 o ] 


g 


13 


= l/3yS 


(2L + 1) (2L 


L ' 

>(2L-1) 2 _ 


1/2 


R 2i t , h 13 - R 10 


1 x + k : 
■ L-l “ 2 


S31 - g 13 [*21~ ^*lo] * ^31 = k 13 [*10 ^*2l] 


, 14 .-/S/rf__kii -l’" R„r‘«,h„ = 

[(2L+1)(2L + 3) 2 J 


*'10 

.L+l 


s 4 1 


“ Sl4 [*2l“ **lo] > k 41 ~ k 14 [*10 _> '*2l] 


S(L, 1 ) 
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r t i 1/2 

g =V3/3 R 21 r L 

(2L+1)(2L-1) 2 _ 


' ^23 ~ ^20 


- + k 2 

1 4 2 s n n 


632 ~ g 23 [^21~ ^^2o] ’ ^32 ~ ^23 [^20 ^ ^ ® 2 l] 


g 2 4=-^3 /3 


1/2 


iii =1 8„ r"',h,.= 


l_(2Ltl)(2L + 3)’J 


20 


2 < r L+l 


§42 ” g 24 [®21*"*®2o] , ^42 “ ^24 [®20~*^2l] 


g 3 4 = - 3 / 3 


L (L + 1) ~ 
[_( 2 L + 1) 4 _ 


1/2 


P r^ + l U _ 21 

k 21 r . n 34 - 

r L 


643 - 634 


In F matrix the interchange of the functions R 10 , R 20 and R 21 
accompanies the interchange of their arguments too. 
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APPENDIX III 


Singularity of the Determinant of Eq. (3.6) 
for L = 0 and 1 

(1) L = 0 Case. By making use of the definition of D Xj and F ij and 
Eq. (3.7) the following relationship can be derived from Eq. (3.3) 


ft 


R 2°l^7 +k H v i - /3rRio ('£i + k *l 


dr 


y^- |®13 ®24 - P 8 2 3 B l+j 


( A5) 


where 


‘13 


= J R 10 R 21 r 3 dr=[2 15 x3- 9 ] 5 


1 2 3 


= j R 20 R 21 r 3 dr = - 3V3 


Integrating the left-hand side of Eq. (A5) by parts, and making use of 
Eqs. (2. 11) and (3. 7) we obtain 


fh«£ 


+ k JjV, “^ rR 10hj^ + k 2j V 2 


dr 


= -2[B„ 


(A6) 


We conclude 


B 21 “ ft B 12 -y^"[ a i3 B 24 “ & 3 2 3 B m] 


(A7) 
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Therefore all the elements on the right-hand side of Eqs. (3.6) are not 
independent. In order that these equations be consistent one of them 
should be a linear combination of others. By making use of the first of 
Eqs. (3.7) it can be shown directly that equations similar to Eq. (A7) 
hold among the elements of each row of the determinant of Eqs. (3.6). 

(11) L = 1 Case. Similar to the previous case the following rela- 
tionship can be derived from Eqs. (3.3) 



where 



Integrating the left-hand side of Eq. (A8) by parts, and making use of 
Eqs. (2. 11) and (3. 7) we obtain 



= -2[B 3 j -/3B 13 ]. ( A9) 
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Combining Eqs. (A8) and (A9) we get 


B 31 -/?B 13 - j^ a i3Bll + /3 a 23Bl2 a 13 ^33 + ^ ^®13 ®34 ~ g^ a 33®14^ 
Finally, Eqs. (3.3) give the following relationship 



Integration by parts of the left-hand side as before gives 



(A10) 


(All) 


= - 2[B 32 - h B 23 ] , (A12) 


whereupon we get 


B 32 - /3 B 23 - — j/3a 23 B 22 + /3a 13 B 21 -a 23 B 33 

+ H ^a 23 B 34 - — /3a 33 B 24 ^ 


(A13) 


To remove the singularity in L = 0 case one of the is chosen 
arbitrary, and a degenerate equation is removed from Eqs. (3.6). Simi- 
larly in L = 1 case two of the are chosen arbitrary and two degenerate 
equations are removed from Eqs. (3.6). 
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APPENDIX IV 


Elements of the Matrix of the Sum of the Asymptotic Coulomb and 

Centrifugal Potentials 

U„ = L(L + 1 ) r' 2 , U 22 = L(L + l)r* 2 

U 33 = (L - 1) L r' 2 + 12(L-1)(2L+1)-V 3 . U 44 = (L+l) (L+2) r~2 + 12 (L+2) (2L+l)- 1 r - 3 
U 12 =U 21 =0, U 13 = U 31 = [256^/243] [L/(2L+1)]V 2 

U, 4 = u 4 ! =- [256T5/243] [(L + l)/( 2L + l)f r’ 3 , U 23 = U 32 = - 6 [l/( 2L + 1)] V 2 
u 24 =U 42 = 6 [(L+1)/(2L+1)] V 2 , U 34 = U 43 =-36 [L(L + 1)] 1/2 (2L + 1) -1 r~ 3 
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TABLES OF CROSS SECTIONS 


In the following five tables the partial and the total cross sections 
for the processes Is-Zs, Is-2p, Is-ls, 2s -2s and 2p-2p in different 
approximations are listed. Each partial cross section corresponds to 
a given total orbital angular momentum L which appears at the head 
of each column. £5 is the sum of the singlet and the sum of the 

triplet partial cross sections within certain approximation and for a 
given electron energy. Q,p is the total cross section, and in the case of 
the eigenstates coupling, Born approximation is used to add the contri- 
bution of the partial waves beyond the maximum L which has been 
computed within the given eigenstates coupling approximation. This 
contribution can easily be obtained by using the given table of the Born 
approximation of the partial and the total cross sections, kj is the wave 
number of the electron beam incident on the ground state of the atom and 
kz is the wave number of the same beam when atom is in its first excited 
state. The energy of the beam is given in electron volt through 

E = 13.6 k2, 

where k can be kj or k 2 . In Table V, lj and l z are the quantum numbers 

of the orbital angular momenta of the bound and the free electrons 
respectively. 
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TABLE I. ls-2s Excitation Cross Sections 


k l 





A. 

Born Approximation 





o 

II 

L = 1 

L = 2 

L = 3 

L =4 

L = 5 

L = 6 

L = 7 

L 

q t 

0.9 

0.16376 


0.00981 

0.00021 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.17378 

0.17379 

1.0 

0.19578 

0.04795 

0.00428 

0.00026 

0.00000 

0.00000 

0.00000 

0.00000 

0.24827 

0.24827 

1.1 

0.16272 

0.07073 

0.01141 

0.00125 

0.00011 

0.00001 

0.00000 

0.00000 

0.24622 

0.24623 

1.2 

0.12704 

0.07896 

0.01858 

0.00299 

0.00039 

0.00004 

0.00000 

0.00000 

0.22800 

0.22800 

1.5 

0.05872 

0.06606 

0.02979 

0.00939 

0.00242 

0.00053 

0.00010 

0.00002 

0.16703 

0.16706 

2.0 

0.01946 

0.03363 

0.02521 

0.01365 

0.00614 

0.00236 

0.00081 

0,00025 

0.10151 

0.10187 

3.0 

0.00388 

0.00909 

0.01019 

0.00866 

0.00628 

0.00394 

0.00226 

0.00116 

0.04546 

0.04758 

4.0 

0.00123 


0.00320 

0.00420 

0.00431 

0.00385 

0.00301 

0.00216 

0.00141 

0.02337 

0.02720 




B. Exchange Neglected 1 

s-2s-2p Eigenstates Coupling Approximation 



o 

II 

>4 

L = 1 

L = 2 

L =3 

n 

L =5 

xO 

II 

►4 

L =7 

L 

Q<p 

0.9 

0.2202 


0.0749 

0.3535 






0.6486 

0.6486 

1.0 

0.1685 


0.1427 

0.1598 

0.0517 





0.5227 

0.5227 

1.1 

0.0951 


0.1142 

0.0298 

0.0616 

0.0231 




0.3238 

0.3238 

1.2 

0.0594 


0.1137 

0.0032 

0.0360 

0.0244 

0.0135 



0.2502 

0.2502 

1.5 

0.0249 


0.0861 

0.0201 

0.0068 

0.0112 

0.0118 

0.0074 


0.1683 

0.1683 

2.0 

0.0101 


0.0373 

0.0255 

0.0107 

0.0046 

0.0034 

0.0033 


0.0949 

0.0953 


C. ls-2s Eigenstates Coupling Approximation 


Singlet 



L = 0 

L = 1 

«M 

II 

L =3 

t- 

ii 

4^ 

L = 5 

L =6 

L =7 

Ss 



0.9 

0.0375 

0.0017 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 


0.0392 



1.0 

0.0725 

0.0583 

0.0002 

0.0000 

0.0000 

0.0000 

0.0000 


0.1310 



1.1 

0.0701 

0.0525 

0.0023 

0.0000 

0.0000 

0.0000 

0.0000 


0.1249 



1.2 

0.0547 

0.0534 

0.0054 

0.0002 

0.0000 

0.0000 

0.0000 


0.1137 



1.5 

0.0241 

0.0384 

0.0110 

0.0022 

0.0004 

0.0001 

0.0000 


0,0762 



2.0 

0.0072 

0.0157 

0.0093 

0.0041 

0.0015 

0.0005 

0.0002 


0.0385 




Triplet 


Qt 












L =0 

L = 1 

L =2 

L =3 

L =4 

L = 5 

L =6 

L =7 




0.9 

0.0004 

0.1686 

0.0060 

0.0000 

0.0000 

0.0000 

0.0000 


0.1750 

0.2142 

0.2142 

1.0 

0.0021 

0.1528 

0.0446 

0.0021 

0.0001 

0.0000 

0.0000 


0.2017 

0.3327 

0.3327 

1.1 

0.0044 

0.1052 

0.0568 

0.0068 

0.0005 

0.0000 

0.0000 


0.1737 

0.2987 

0.2987 

1.2 

0.0061 

0.0737 

0.0576 

0.0114 

0.0015 

0.0002 

0.0000 


0.1505 

0.2642 

0.2642 

1.5 

0.0073 

0.0355 

0.0406 

0.0174 

0.0050 

0.0012 

0.0002 


0.1072 

0.1833 

0.1833 

2.0 

0.0049 

0.0162 

0.0205 

0.0143 

0.0074 

0.0032 

0.0012 


0.0677 

0.1062 

0.1068 




D. Is- 

2s -2p Eigenstates 

Coupling Approximation 





Singlet 




L = 0 

L = 1 

L = 2 

L = 3 

L =4 

L = 5 

L =6 

L =7 




0.9 

0.0523 

0.0048 

0.0620 






0.1191 



1.0 

0.0768 

0.0147 

0.0833 

0.0092 





0.1840 



1.1 

0.0585 

0.0245 

0.0647 

0.0236 

0.0055 




0.1768 



1.2 

0.0382 

0.0251 

0.0246 

0.0252 

0.0081 

0.0028 



0,1240 



1.5 

0.0123 

0.0308 

0.0015 

0.0041 

0.0051 

0.0034 

0.0023 

0.0026 

0.0621 



2.0 

0.0049 

0.0152 

0.0068 

0.0021 

0.0010 

0.0010 

0.0011 

0.0008 

0.0329 



3.0 

0.0010 

0.0031 

0.0031 

0.0023 

0.0015 

0.0009 

0.0005 

0.0006 

0.0130 



4.0 

0.0003 

0.0010 

0.0012 

0.0012 

0.0010 

0.0008 



0.0055 




Triplet 

2$+2 t 

Qj 












L =0 

L = 1 

L = 2 

L = 3 

L = 4 

II 

L = 6 

L = 7 

2 T 



0.9 

0.0013 

0.0748 

0.0019 






0.0780 

0.1971 

0.1971 

1.0 

0.0040 

0.1224 

0.0195 

0.0214 





0.1673 

0.3513 

0.3513 

1.1 

0.0050 

0.1013 

0.0326 

0.0077 

0.0131 




0.1597 

0.3366 

0.3366 

1.2 

0.0055 

0.0724 

0.0359 

0.0036 

0.0105 

0.0076 



0.1355 

0.2596 

0.2596 

1.5 

0.0045 

0.0333 

0.0309 

0.0072 

0.0046 

0.0054 

0.0049 

0.0043 

0.0951 

0.1573 

0.1573 

2.0 

0.0031 

0.0155 

0.0176 

0.0101 

0.0044 

0.0025 

0.0023 

0.0019 

0.0574 

0.0903 

0.0907 

3.0 

0.0013 

0.0048 

0.0065 

0.0059 

0.0044 

0.0029 

0.0018 

0.0012 

0.0288 

0.0418 

0.0439 

4.0 

0.0006 

0.0019 

0.0028 

0.0030 

0.0027 

0.0022 



0.0132 

0.0187 

0.0261 


54 


TABLE II. ls-2p Excitation Cross Sections 



0.00107 

0.00499 

0.00702 

0.00747 

0.00550 

0.00234 

0.00048 

0.00013 


0.46700 | 
0.48867 ' 
0.36207 ! 
0.25540 
0.09184 
0.02222 
0.00285 
0.00066 


0.09607 

0.35645 

0.41698 

0.38137 

0.19986 

0.05938 

0.00806 

0.00179 


0.01025 

0.13467 

0.24992 

0.29903 

0.23775 

0.09386 

0.01537 

0.00359 


0.00087 

0.03964 

0.11533 

0.17908 

0.21386 

0.11287 

0.02309 

0.00586 


0.00007 1 
0.01032 
0.04649 
0.09269 
0.16292 
0.11376 
0.02920 
0.00816 


0.00001 

0.00254 

0.01753 

0.04453 

0.11291 

0.10248 

0.03271 

0.01009 


0.00000 

0.00059 

0.00624 

0.02004 

0.07232 

0.08374 

0,03273 

0.01114 


0.56534 

1.03787 

1.22158 

1.27961 

1.09696 

0.59065 

0.14449 

0.04142 


0.57535 

1.03851 

1.22859 

1.30741 

1.28101 

1.04055 

0.66256 

0.45252 



B. Exchange Neglected ls-2s-2p Eigenstates Coupling Approxir 
























TABLE III. Is- Is Elastic Cross Sections 






A. 

Born Approximation 





o 

II 

L = 1 

L = 2 

L s 3 

L = 4 

L = 5 

L = 6 

it 

L 

Q <2* 

0.9 

1.65113 

0.09722 

0.00379 

0.00013 

0.00000 

0.00000 

t 

0.00000 

0.00000 

1.75227 

1.75233 

1.0 

1.42354 

0.11192 

0.00587 

0.00027 

0.00001 

0.00000 

0.00000 

0.00000 

1.54161 

1.54167 

1.1 

1.22728 

0.12348 

0.00835 

0.00049 

0.00003 

0.00000 

0.00000 

0.00000 

1.35963 

1.35968 

1.2 

1.05936 

0.13181 

0.01109 

0.00082 

0.00006 

0.00000 

0.00000 

0.00000 

1.20314 

1.20318 

1.5 

0.69140 

0.14017 

0.01955 

0.00241 

0.00028 

0.00003 

0.00000 

0.00000 

0.85384 

0.85389 

2.0 

0.36283 

0.12285 

0.02924 

0.00621 

0.00124 

0.00023 

0.00004 

0.00001 

0.52265 

0.52267 

3.0 

0.12662 

0.07286 

0.03013 

0.01131 

0.00402 

0.00133 

0.00042 

0.00012 

0.24681 

0.24700 

4.0 

0.05565 

0.04195 

0.02291 

0.01145 

0.00546 

0.00243 

0.00104 

0.00042 

0.14131 

0.14194 



B. Exchange Neglected Is- 

2s -2p Eigenstates Coupling Approximation 



L= 0 

L = 1 

L = 2 

L = 3 

L = 4 

L= 5 

L = 6 

L = 7 

I 

Qx 

0.9 

2.632 

0.3300 

0.1916 

0.0157 





3.169 

3.169 

1.0 

2.091 

0.2518 

0.0890 

0.0309 

0.0078 

0.0028 



2.473 

2.473 

1.1 

1.694 

0.2302 

0.0483 

0.0299 

0.0105 

0.0038 

0.0016 


2.018 

2.018 

1.2 

1.390 

0.2228 

0.0346 

0.0248 

0.0113 

0.0045 

0.0019 


1.690 

1.690 

1.5 

0.8181 

0.2032 

0.0310 

0.0135 

0.0091 

0.0048 

0.0027 


1.082 

1.082 

2.0 

0.4053 

0.1742 

0.0494 

0.0181 

0.0108 

0.0036 

0.0027 


0.6641 

0.6641 


C. Is -2s Eigenstates Coupling Approximation 


Singlet 



o 

ii 

-3 

L =1 

cq 

II 

L =3 

n 

4 ^ 

r 

II 

U 9 

L =6 

L = 7 




0.9 

0.4324 

0.0382 

0.0015 

0.0000 

0.0000 

0.0000 

0.0000 


0.4721 



1.0 

0.2824 

0.0338 

0.0014 

0.0001 

0.0000 

0.0000 

0.0000 


0.3177 



1.1 

0.1865 

0.0199 

0.0011 

0.0001 

0.0000 

0.0000 

0.0000 


0.2076 



1.2 

0.1397 

0.0111 

0.0006 

0.0001 

0.0000 

0.0000 

0.0000 


0.1515 



1.5 

0.0905 

0.0012 

0.0001 


0.0000 

0.0000 

0.0000 


0.0918 



2.0 

0.0608 

0.0083 

0.0016 

0.0003 

0.0001 

0.0000 

0.0000 


0.0711 




Triplet 

2$ +2 t 









— 



Qt 


r 

n 

o 

L =1 

t- 

ii 

t\> 

L = 3 

II 

a 

in 

II 

L =6 

L= 7 




0.9 

3.684 

1.394 

0.0477 

0.0016 

0.0001 

0.0000 

0.0000 


5.127 

5.599 

5.599 

1.0 

2.903 

1.162 

0.0579 

0.0023 

0.0001 

0.0000 

0.0001 


4.125 

4.443 

4.443 

1.1 

2.297 

0.9654 

0.0651 

0.0036 

0.0002 

0.0000 

0.0000 


3.331 

3.539 

3.539 

1.2 

1 .829 

0.8085 

0.0696 

0.0050 

0.0004 

0.0000 

0.0000 


2.713 

2.864 

2.864 

1.5 

0.9716 

0.4857 

0.0727 

0.0092 

0.0011 

0.0001 

0.0000 


1.5404 

1.632 

1.632 

2.0 

0.4042 

0.2366 

0.0627 

0.0142 

0.0031 

0.0006 

0.0001 


0.7215 

0.7926 

0.7926 
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TABLE III. Continued 


k l 



D. 

ls-2p Eigenstates Coupling Approximation 




Singlet 




L =0 

L =1 

L = 2 

L = 3 

L = 4 

L = 5 

L = 6 

L = 7 

VI 



0.9 

0.4230 

0.0186 

0.0442 

0.0043 





0.4901 



1.0 

0.2984 

0.0167 

0.0635 

0.0073 

0.0017 

0.0006 

0.0005 


0.3887 



1.1 

0.2148 

0.0108 

0.0578 

0.0104 

0.0026 

0.0008 

0.0003 


0.2975 



1.2 

0.1569 

0.0048 

0.0338 

0.0105 

0.0030 

0.0011 

0.0005 


0.2106 



1.5 

0.1040 

0.0008 

0.0035 

0.0038 

0.0025 

0.0012 

0.0007 


0.1165 



2.0 

0.0698 

0.0138 

0.0030 

0.0018 

0.0018 

0.0008 

0.0006 


0.0916 




Triplet 

2$ +2 t 

Qt 


L = 0 

L= 1 

L= 2 

L = 3 

L = 4 

L=5 

L = 6 

L = 7 


0.9 

3.689 

1.421 

0.0682 

0.0136 





5.192 

5.682 

5.682 

1.0 

2.914 

1.236 

0.0758 

0.0152 

0.0054 

0.0019 

0.0015 


4.250 

4.638 

4.638 

1.1 

2.309 

1.036 

0.0814 

0.0142 

0.0066 

0.0025 

0.0011 


3.451 

3.749 

3.749 

1.2 

1.840 

0.8702 

0.0850 

0.0139 

0.0069 

0.0031 

0.0015 


2.821 

3.031 

3.031 

1.5 

0.9806 

0.5249 

0.0863 

0.0152 

0.0063 

0.0031 

0.0019 


1.618 

1.735 

1.735 

2.0 

0.4184 

0.2702 

0.0828 

0.0256 

0.0111 

0.0030 

0.0019 


0.8130 

0.9045 

0.9045 


E. ls-2s-2p Eigenstates Coupling Approximation 




Singlet 




L =0 

L =1 

L =2 

L =3 

L =4 

L= 5 

L =6 

L= 7 

2 S 



0.9 

0.4503 

0.0133 

0.0457 

0.0058 

0.0000 




0.5151 



1.0 

0.2649 

0.0100 

0.0660 

0.0076 

0.0019 

0.0007 

0.0007 


0.3518 



1.1 

0.1736 

0.0131 

0.0582 

0.0104 

0.0028 

0.0009 

0.0004 


0.2594 



1.2 

0.1283 

0.0102 

0.0351 

0.0104 

0.0030 

0.0012 

0.0005 


0.1887 



1.5 

0.0854 

0,0021 

0.0049 

0.0042 

0.0025 

0.0012 

0.0007 

0.0010 

0.1020 



2.0 

0.0591 

0.0093 

0.0016 

0.0009 

0.0007 

0.0006 

0.0004 

0.0002 

0.0728 



3.0 

0.0266 

0.0138 

0.0052 

0.0019 

0.0008 

0.0003 

0.0002 

0.0005 

0.0493 



4.0 

0.0133 

0.0106 

0.0061 

0.0034 

0.0020 

0.0012 



0.0366 




Triplet 

2$ +2 t 

Q ■p 


n 

o 

L = 1 

L= 2 

L =3 

L =4 

L =5 

L =6 

L =7 

S T 

0.9 

3.687 

1.726 

0.0875 

0.0166 

0.0055 




5.523 

6.037 

6.037 

1.0 

2.908 

1.372 

0.0938 

0.0174 

0.0062 

0.0021 

0.0016 


4.401 

4.753 

4.753 

1.1 

2.299 

1.090 

0.0949 

0.0166 

0.0075 

0.0028 

0.0012 


3.512 

3.772 

3.772 

1.2 

1.828 

0.8884 

0.0956 

0.0162 

0.0076 

0.0033 

0.0015 


2.841 

3.029 

3.029 

1.5 

0.9697 

0.5228 

0.0911 

0.0166 

0.0069 

0.0033 

0.0021 

0.0020 

1.615 

1.716 

1.716 

2.0 

0.4006 

0.2433 

0.0657 

0.0154 

0.0044 

0.0020 

0.0011 

0.0008 

0.7333 

0.8061 

0.8061 

3.0 

0.1127 

0.0854 

0.0388 

0.0156 

0.0062 

0.0022 

0.0009 

0.0007 

0.2625 

0.3118 

0.3120 

4.0 

0.0475 

0.0444 

0.0276 

0.0160 

0.0093 

0.0053 



0.1501 

0.1867 

0.1888 
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Born Approximation 



L = 0 

L = 1 

L = 2 

L = 3 

L = 4 

L = 5 

L = 6 

L = 7 

I 

Qt 

0.24 

389.97 

15.910 

0.27680 

0.00319 

0.00003 

0.00000 

0.00000 

0.00000 

406.16 

406.17 

0.50 

105.52 

32.129 

5.0152 

0.55306 

0.04846 

0.00347 

0.00022 

0.00001 

143.269 

143.276 

0.68 

47.381 

24.746 

7.4981 

1.7023 

0.31602 

0.04869 

0.00658 

0.00078 

81.700 

81.703 

0.83 

26.562 

17.895 

7.5601 

2.5080 

0.69820 

0.16365 

0.03393 

0.00621 

55.427 

55.440 

1.23 

8.2461 

7.6846 

4.9702 

2.7483 

1.3500 

0.57745 

0.22452 

0.07803 

25.879 

25.990 

1.80 

2.4331 

2.7773 

2.2763 

1.6805 

1.1541 

0.71994 

0.41776 

0.22120 

11.680 

12.105 

2.87 

0.5245 

0.71381 

0.68903 

0.61123 

0.51876 

0.41048 

0.31076 

0.21860 

4.0011 

4.8280 

3.91 

0.18794 

0.27611 

0.28658 

0.27339 

0.25098 

0.21580 

0.17946 

0.14029 

1.8106 

2.7417 


B. Exchange Neglected ls-2s-2p Eigenstates Coupling Approximation 



L ^ 0 

L - 1 

L = 2 

L = 3 

L = 4 

L = 5 

L = 6 

L =7 £ 

Q x 

0.24 

32.02 

8.489 

42.26 





82.77 

82.78 

0.50 

2.041 

7.710 

23.35 

9.004 

9.368 

7.702 

5.870 

65.045 

65.052 

0.68 

2.255 

8.247 

15.16 

5.547 

4.373 

3.755 

2.975 

42.312 

42.316 

0.83 

1.861 

8.206 

10.67 

4.637 

2.792 

2.223 

1.818 

32.207 

32.226 

1.23 

1.716 

4.987 

5.030 

3.321 

1.797 

1.058 

0.7179 

18.537 

18.726 

1.80 

1.020 

2.134 

2.150 

1.750 

1.275 

0.8320 

0.5457 

9.707 

10.352 


C. ls-2s Eigenstates Coupling Approximation 
Singlet 



L = 0 

L - 1 

L = 2 

L= 3 

L = 4 

L = 5 

L = 6 

L = 7 2 S 

0.24 

0.3303 

8.196 

0.2628 

0.0028 




8.792 

0.50 

1.532 

10.38 

0.0275 

0.0048 

0.0008 

0.0002 

0.0249 

11.97 

0.68 

1.115 

5.536 

1.502 

0.1150 

0.0087 

0.0010 

0.0017 

8.279 

0.83 

0.8980 

3.512 

1.997 

0.4303 

0.0747 

0.0129 

0.0032 

6.928 

1.23 

0.5702 

1.413 

1.236 

0.7010 

0.3129 

0.1228 

0.0450 

4.401 

1.80 

0.2825 

0.5370 

0.5285 

0.4193 

0.2931 

0.1863 

0.1110 

2.358 



— 



Triple 





v + v 


L - 0 

L - 1 

L = 2 

L = 3 

L = 4 

L = 5 

L = 6 

L = 7 

2t 


Q x 

45.94 

0.2102 

118.8 

34.44 

7.713 

21.05 

0.0540 

2.776 

0.2521 

0.0316 

0.0994 


172.51 

58.86 

181.30 

70.83 

181.31 

70.84 

1.366 

18.13 

12.74 

4.059 

0.8282 

0.1463 

0.0442 


37.31 

45.60 

45.60 

2.112 

11.65 

8.725 

3.887 

1.230 

0.3225 

0.0850 


28.01 

34.94 

34.96 

1.811 

4.691 

4.008 

2.585 

1.399 

0.6637 

0.2862 


15.44 

19.84 

20.03 

0.8989 

1.735 

1.652 

1.316 

0.9510 

0.6359 

0.4008 


7.590 

9.947 

10.592 


D. ls-2s-2p Eigenstates Coupling Approximation 
Singlet 



L = 0 

L. = 1 

L - 2 

L = 3 

L - 4 

L = 5 

0.24 

7.800 

14.79 

22.42 




0.50 

0.2858 

0.6960 

3.149 

4.447 

2.925 

2.063 

0.68 

0.0661 

1.044 

2.455 

1.884 

1.480 

1.098 

0.83 

0.1675 

1.088 

2.105 

0.9905 

0.8282 

0.6606 

1.23 

0.3739 

1.075 

1.196 

0.7196 

0.3940 

0.2496 

1.80 

0.2416 

0.4974 

0.5133 

0.4133 

0.2923 

0.1928 

2.87 

0.0847 

0.1489 

0.1574 

0.1466 

0.1280 

0.1057 

3.91 

0.0365 

0.0635 

0.0702 

0.0698 

0.0660 

0.0599 


Triplet 


L = 0 

L. = 1 

L - 2 

L - 3 

L = 4 

L = 5 


0.24 I6.52 a 

0.50 6.172 


1.236 

17.19 

12.88 

9.166 

4.199 

1.656 

0.4652 

0.1945 


89.65 

30.69 

14.52 

8.976 

3.880 

1.603 

0.4857 

0.2146 


12.90 

7.937 

5.367 

2.681 

1.291 

0.4479 

0.2122 


45.01 

15.057 

8.820 

6.347 

4.332 

2.3553 

0.9223 

0.3659 


107.41 

82.14 

45.28 

30.343 

15.622 

7.648 

2.8413 

1.1138 


152.42 

97.206 

54.099 

36.686 

19.955 

10.002 

3.7636 

1.4797 


TABLE V. 2p-2p Elastic Cross Sections 


k 2 


A. 

L-V*2 

Odd, Born Approximation 





L = 2 

L = 3 

L = 4 

L = 5 

L = 6 

L = 7 

2 0 

0.24 

26.562 

6.1961 


0.92680 

0.44497 

0.22694 

0.08105 

36.56 

0.50 

14.210 

5.2190 

2.1768 

1.0526 

0.55433 

0.31476 

0.18383 

23.71 

0.68 

8.8346 

4.0761 

1.9476 

1.0100 

0.54983 

0.31728 

0.18718 

16.923 

0.83 

6.0065 

3.2062 

1.6993 

0.94122 

0.53204 

0.31331 

0.18680 

12.8854 

1.23 

2.5053 

1.7048 

1.1013 

0.70920 

0.44518 

0.28212 

0.17577 

6.9237 

1.80 

0.91169 

0.74615 

0.56949 

0.42475 

0.30158 

0.21122 

0.14245 

3.30733 

2.87 

0.23647 

0.22647 

0.20035 

0.17178 

0.13878 

0.10926 

0.08177 

1.16488 


B. L-4, - 

Ji 2 Odd, Exchange Neglected 2] 

p Eigenstates Coupling Approximations 


■SB 

L = 2 

— ■ 

m 

L = 5 

■m 

F 

ii 

-4 

2 0 

0.24 

61.12 

8.444 

2.5808 

1.1408 

0.6200 

0.3720 

0.2516 

74.52 

0.50 

15.292 

6.884 

2.5436 

1.1580 

0.6216 

0.3640 

0.2352 

27.10 

0.68 

8.008 

4.940 

2.2548 

1.1160 

0.6160 

0.3640 

0.2352 

17.54 

0.83 

5.108 

3.6364 

1.9292 

1.0380 

0.5968 

0.3592 

0.2340 

12.90 

1.23 

2.0484 

1.7392 

1.1812 

0.7652 

0.4972 

0.3236 

0.2208 

6.776 

1.80 

0.7640 

0.7220 

0.5816 

0.4436 

0.3312 

0.2400 

0.1784 

3.261 

2.87 

0.2132 

0.2180 

0.2004 

0.1776 

0.1520 

0.1272 

0.1064 

1.195 



C. L-4, 

-4 2 Odd, Zp Eigenstates Coupling Approximation 



Singlet 


L = 1 

L = 2 

L = 3 

L = 4 

L = 5 

L = 6 

L = 7 

2 os 

0.24 

2.963 

4.161 

0.6725 

0.2861 

0.1552 

0.0930 

0.0629 

8.394 

0.50 

3.735 

3.182 

0.7915 

0.3066 

0.1576 

0.0913 

0.0589 

8.323 

0.68 

2.165 

1.728 

0.6851 

0.3025 

0.1587 

0.0920 

0.0590 

5.190 

0.83 

1.371 

1.107 

0.5562 

0.2799 

0.1547 

0.0913 

0.0589 

3.619 

1.23 

0.5280 

0.4603 

0.3115 

0.1991 

0.1277 

0.0823 

0.0558 

1.7647 

1.80 

0.1928 

0.1831 

0.1475 

0.1124 

0.0837 

0.0605 

0.0449 

0.8249 

2.87 

0.0534 

0.0547 

0.0503 

0.0445 

0.0381 

0.0319 

0.0267 

0.2996 


Triplet 


L = 1 

L = 2 

L = 3 

L = 4 

L = 5 

L = 6 

L = 7 

S 0T 

0.24 

49.22 

3.265 

1.850 

0.8528 

0.4650 

0.2788 

0.1886 

56.12 

0.50 

7.791 

1.851 

1.504 

0.8193 

0.4598 

0.2722 

0.1764 

12.87 

0.68 

4.720 

2.219 

1.361 

0.7700 

0.4481 

0.2700 

0.1756 

9.964 

0.83 

3.373 

2.099 

1.237 

0.7198 

0.4318 

0.2652 

0.1744 

8.300 

1.23 

1.484 

1.225 

0.8376 

0.5509 

0.3630 

0.2384 

0.1639 

4.863 

1.80 

0.5674 

0.5337 

0.4295 

0.3283 

0.2456 

0.1784 

0.1330 

2.4159 

2.87 

0.1596 

0.1632 

0.1500 

0.1328 

0.1136 

0.0951 

0.0797 

0.8940 
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TABLE V. (Continued) 


k 2 

D. L Even, Born Approximation 

0.24 

0.50 

0.68 

0.83 

1.23 

1.80 

2.87 

3.91 

0.24 

0.50 

0.68 

0.83 

1.23 

1.80 

0.24 

0.50 

0.68 

0.83 

1.23 

1.80 

0.24 

0.50 

0.68 

0.83 

1.23 

1.80 

0.24 

0.50 

0.68 

0.83 

1.23 

1.80 

2.87 

3.91 

0.24 

0.50 

0.68 

0.83 

1.23 

1.80 

2.87 

3.91 

L = 0 

L = 1 

L = 2 

L = 3 

L = 4 

L = 5 

L = 6 

L = 7 

S E 

Qt 

12.488 

0.1758 

0.07386 

0.22032 

0.25024 

0.13477 

0.04395 

0.01868 

230.42 

74.475 

36.478 

21.559 

7.3349 

2.3518 

0.55809 

0.20876 

4.0427 

7.0963 

6.5573 

5.4143 

2.9703 

1.3155 

0.40443 

0.17048 

1.9014 

0.92172 

1.0156 

1.1653 

1.1453 

0.74676 

0.30267 

0.14307 

0.30797 

0.48094 

0.29902 

0.26124 

0.37750 

0.39290 

0.21866 

0.11656 

0.57480 

0.45745 

0.27895 

0.15835 

0.11091 

0.18754 

0.15188 

0.09255 

0.29675 

0.30098 

0.23113 

0.15225 

0.04377 

0.07731 

0.09753 

0.06896 

0.06128 

0.05695 

0.04364 

0.02533 

0.03572 

0.02824 

0.05832 

0.04864 

250.09 

83.965 

44.978 

28.956 

12.269 

5.2348 

1.8355 

0.8677 

290.69 

109.714 

63.742 

43.572 

20.571 

9.6066 

3.8272 

2.1455 

E. L-'t l -'t 2 Even, Exchange Neglected ls-2s-2p Eigenstates Coupling Approximation 

L = 0 

L = 1 

L = 2 

L = 3 

L = 4 

L = 5 

L = 6 

L = 7 


Qt 

31.96 

9.371 

4.156 

2.542 

1.208 

.5612 

91.21 

12.65 

5.953 

4.560 

2.879 

1.385 

154,6 

30.88 

13.19 

7.424 

2.906 

1.238 

92.86 

20.53 

9.628 

5.249 

1.734 

.8299 

55.86 

12.99 

6.563 

3.811 

1.106 

.5109 

36.79 

8.451 

4,508 

2.793 

.8438 

.3005 

24.37 

5.695 

3.130 

2.032 

.6912 

.2022 


487.7 

100.6 

47.13 

28.41 

11.37 

5.028 

566.32 

129.8 

66.56 
43.07 

19.56 
9.382 

F. L-^j -^2 Even, ls-2p Eigenstates Coupling Approximation 


Singlet 


L = 0 

L = 1 

L = 2 

L = 3 

L = 4 

L = 5 

L = 6 

L = 7 

Z ES 

1.964 

0.5131 

0.2346 

0.1227 

0.1071 

0.0571 

5.238 

3.159 

2.091 

1.252 

0.6796 

0.3366 

17.34 

1.544 

1.039 

0.9930 

0.6167 

0.2926 

1.260 

1.439 

0.4944 

0.1532 

0.2249 

0.1818 

0.4896 

0.4063 

0.2539 

0.1191 

0.0635 

0.0954 

0.2460 

0.1906 

0.1520 

0.0990 

0.0249 

0.0426 

0.1133 

0.0984 

0.0766 

0.0217 

0.0172 


26.54 

7.365 

4.363 

2.816 

1.7384 

1.0233 

Triplet 

Qt 

L = 0 

L = 1 

L = 2 

L = 3 

L = 4 

L = 5 

L = 6 

L = 7 

2 et 

15.62 

4,844 

2.247 

1.293 

0.4842 

0.1909 

26.72 

3.384 

3.886 

3.650 

2.287 

1.063 

32.44 

19.97 

10.33 

6.287 

2.424 

0.9673 

6.674 

12.52 

5.922 

3.476 

1.418 

0.6659 

1.252 

0.8285 

1.035 

0.9248 

0.6072 

0.4001 

0.7359 

0.3266 

0.2270 

0.2097 

0.2079 

0.2024 

0.2754 

0.1764 

0.1162 

0.0717 

0.0962 


83.44 

42.15 

23,82 

15,957 

7,500 

3.586 

178.89 

72.80 

45.22 

32.45 

17.28 

8.943 

G. Even, ls-2s-2p Eigenstates Coupling Approximation 


Singlet 

L = 0 

L = 1 

L = 2 

L = 3 

L = 4 

L = 5 

L = 6 

L = 7 


7.852 

2.470 

1.344 

.7424 

.2813 

.1357 

.0513 

.0267 

13.45 

5.026 

2.283 

1.316 

.6752 

.3339 

.1053 

.0451 

38.56 

7.433 

3.025 

1.579 

0.6254 

.2862 

.0930 

.0422 

21.41 

4.900 

2.683 

1.518 

.3663 

.1816 

.0739 

.0373 

15.44 

3.201 

1.756 

1.079 

.2726 

.1065 

.0553 

.0318 

8.610 

2.101 

1.166 

.7580 

.2289 

.0659 

.0391 

.0261 

6.432 

1.417 

.7903 

.5313 

.1916 

.0471 

.0268 

.0212 

.0383 

.0181 

.0166 

111.75 

26.548 

13.047 

7.524 

2.641 

1.1952 

0.4628 

0.2470 

Triplet 

L = 0 

L = 1 

L = 2 

L = 3 

L = 4 

L = 5 

L = 6 

L = 7 

S ET 

Q t 

27.90 

4.337 

2.823 

2.032 

.9741 

.4291 

.1564 

.0805 

75.79 

4.540 

4.018 

3.674 

2.290 

1.055 

.3217 

.1366 

63.10 

20.66 

10.75 

6.400 

2.465 

.9303 

.2907 

.1285 

87.21 

21.62 

8.986 

4.796 

1.599 

.6614 

.2330 

.1143 

41.31 

10.38 

5.149 

2.950 

.9568 

.4150 

.1780 

.0982 

26.12 

6.373 

3.337 

2.013 

.6421 

.2540 

.1291 

.0810 

17.13 

4.255 

2.302 

1.468 

.4924 

.1639 

.0907 

.0655 

.1180 

.0623 

.0524 

338.56 

72.17 

37.37 

23.333 

9.419 

4.027 

1.4619 

0.7570 

518.92 

122.01 

67.46 

44.54 

20.11 

9.531 

3.945 
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Figure 1 — ls-2s Excitation Cross Section. The experimental curves 
are given in references 1 and 2. The ls-2s curve corre- 
sponds to a wave function which contains the Is and the 2s 
atomic states. The Is-2s-2p curve corresponds to the ad- 
dition of the 2p states to the above wave function. The 
(3 = 0 curve is the case when the exchange potentials are 
neglected. In ls-2s-2p sym it is assumed that electrons 
are spinless particles whose total wave function is symmetric. 

Figure 2 — ls-2p Excitation Cross Section. The experimental curve is 
given in reference 3. The ls-2p curve corresponds to a 
wave function which contains the Is and the 2p atomic states. 
The rest of the curves have the same meaning as in Figure 1. 

Figure 3 -- Is- Is Elastic Cross Section. The curves have the same 
meaning as in Figure 1 and 2. 

Figure 4 — 2s-2s Elastic Cross Section. The curves have the same 
meaning as in Figures 1 and 2. 

Figure 5 — 2p-2p Elastic Cross Section. The curves have the same 
meaning as in Figure 2. 
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FIG. I. Is -2s EXCITATION CROSS SECTION 
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FIG. 2. ls-2p EXCITATION CROSS SECTION 
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FIG. 3. Is-ls ELASTIC CROSS SECTION 
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FIG. 5. 2p-2p ELASTIC CROSS SECTION 
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